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We consider the relationship between the temperature at which averaged energy landscape proper- 
ties change sharply (T ), and the breakdown of mean-field treatments of the dynamics of supercooled 
liquids. First, we show that the solution of the wavevector dependent mode-coupling equations un- 
dergoes an ergodic-nonergodic transition consistently close to T . Generalizing the landscape concept 
to include hard-sphere systems, we show that the property of inherent structures that changes near 
To is governed more fundamentally by packing and free volume than potential energy. Lastly, we 
study the finite-size Random Orthogonal Model (ROM), and show that the onset of noticeable cor- 
rections to mean-field behavior occurs at T . These results highlight new connections between the 
energy landscape and mode-coupling approach to supercooled liquids, and identify what features of 
the relaxation of supercooled liquids are properly captured by mode-coupling theory. 
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The slow dynamics that supercooled liquids undergo 
as they approach the glass transition has defied a satisfy- 
ing explanation for many decades^. Several theoretical 
paradigms have been presented that shed light on cer- 
tain features of these dynamics. The notion of an energy 
landscape has been useful for understanding thermody- 
namic properties of the glassy state as well as for ratio- 
nalizing the connection between transport properties and 
the packing structures associated with local potential en- 
ergy minima ("inherent structures") that are visited at 
a given temperature^. Other theories, such as mode- 
coupling theory (MCT)2£, have also been influential in 
explaining the sequence of relaxation events that occurs 
in mildly supercooled liquids. While successful in several 
contexts, both the landscape and MCT approaches suf- 
fer from problems that limit their utility. The landscape 
picture does not provide a predictive, quantitative frame- 
work for describing the dynamics of supercooled liquids. 
Furthermore, it is not clear how one could apply energy 
landscape concepts in a useful way to entropically dom- 
inated glassy systems such as hard-sphere liquids. MCT 
does provide such a framework, but several important 
predictions made by this theory, including the thermody- 
namic location of the ergodic-nonergodic transition, are 
incorrect. Here we elucidate connections between these 
two viewpoints that shed new light on various features of 
these seemingly different approaches. 

Effort has been made to connect the intuitively based 
landscape picture to the more mathematical mode- 
coupling approach. Pioneering work of Kirkpatrick, 
Wolynes and Thirumalai^ showed that the mode- 
coupling equations are exact for certain mean-field spin 
glasses. In particular, the p — 3 p-spin model exhibits a 
dynamical transition at a temperature T c and a thermo- 
dynamic transition at a lower temperature 7 1 R J2 4 Ifi i ii i i£. 
The temperature T c is associated with a mean-field di- 
vergence of barriers, leading to nonergodicity. It has 
been argued that effects beyond the mean-field limit ren- 
der the barriers at T c finite*^. Thus, T c corresponds 



to the temperature at which activated processes over fi- 
nite sized barriers dominate transport. Since MCT in its 
simplest form neglects these activated processes, the full 
wavevector dependent MCT equations exhibit a "glass" 
transition at T c > T g where T g is the calorimetric glass 
transition observed in the laboratory. Several computer 
studies have attempted to strengthen the connection be- 
tween activated processes on the energy landscape and 
the temperature Recently, however, simula- 

tions have shown that activated processes strongly in- 
fluence transport in supercooled liquids at temperatures 
significantly in excess of T,mLl^ii^. Thus, the physical 
meaning and relevance of T c for finite dimensional, non- 
mean-field glassy systems remains unclear. 

Several years ago, Sastry et al& pointed out the ex- 
istence of another characteristic landscape temperature, 
the "onset" (or "landscape influenced") temperature T a . 
T Q may be significantly larger than T c as calculated by 
power-law fits of diffusive data, and coincides with the 
onset of nonexponential and non-Arrhenius relaxation in 
supercooled liquids. Sastry et al. 3 found that T also 
marked the temperature at which averaged energy land- 
scape properties (such as the average value of the inher- 
ent structure energy) show a sharp change as a function 
of temperature. While much work has focused on un- 
derstanding the qualitative changes in dynamics near T c , 
very little work has been devoted to understanding the 
meaning of T Q . Here we demonstrate connections be- 
tween MCT, (and mean-field-like approaches in general) 
and the changes in landscape properties that occur as the 
system is cooled below T a . 

It has long been known that the location of the temper- 
ature at which the full wavevector dependent MCT equa- 
tions predict a loss of ergodicity is significantly higher 
than T c as obtained by power-law fits of the temperature 
dependence of transport coefficients. Here we show that 
the ergodic-nonergodic transition actually occurs close to 
T a . We solve the wavevector dependent MCT equations 
for the 2x2 matrix F(q, t) with matrix elements Fij (q, i), 
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where 1 ' 



F(q,i) + n 2 (q)F(q,i)+ / drM(q, t - r)F(q, r) = 0. 



The frequency matrix l~i 2 (q) is defined as 



(1) 



E^[ s_1 (q)]^ ( 2 ) 



and the mode-coupling approximation for the memory 
function is given by 

^•a'/J' (q, q - k)F aQ , (k, tJF/j/j/ (q - k, t) (3) 

where n is the particle density, Xi and rrij are, respec- 
tively, the concentration and mass of particle type i and 
S(q) is the 2x2 matrix of partial structure factors. 
Expressions for the vertices Vi a /3 can be found in the 
literature^. 
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FIG. 1: Average inherent structure energy vs. temperature 
for (a) 70 particle 50/50 soft-sphere system and (b) 120 par- 
ticle 80/20 LJ system. The solid lines are straight line fits 
to high and low temperature slopes. T is defined to be the 
temperature where these lines cross. The vertical dashed line 
indicates the temperature of the ergodic-nonergodic transi- 
tion as calculated by Eqs. (l)-(3). T c , as calculated by fits to 
diffusive data, is found to be T c = 0.246 for the soft-sphere 
system and T c = 0.435 for the LJ system. Lennard- Jones 
units are used here and throughout the text. Both systems 
are simulated at unit density. 

Solving Eqs. (l)-(3) yields the location of the temper- 
ature at which the function F(\a,t)/S(k) fails to decay 
to zero as t — > oo. In Fig. 1 we plot the average in- 
herent structure energy (Eis(T)) as a function of tem- 
perature for two different potentials. The first system 
studied is the 50/50 soft-sphere mixture studied by Bar- 
rat et. alM. The second system is the 80/20 Lennard- 
Jones mixture studied by Kob and AndersenSi. The 
onset temperature T Q is located where Ej${T) show a 
sharp decrease as a function of temperature, or where 
the 1st order polynomial fits to the high and low temper- 
ature data cross. A dashed line indicates the location of 



the ergodic-nonergodic transition temperature obtained 
from Eqs. (l)-(3). Clearly, the location of the ergodic- 
nonergodic transition occurs very close to the landscape 
onset temperature T Q and not T c as calculated by fits to 
diffusion data. Since MCT may be viewed as a particular 
type of dynamical mean-field theory, the coincidence of 
the breakdown of MCT at T D signals the failure of MCT 
to capture specific non-mean-ficld effects. The nature of 
this failure is discussed below. 

A particularly troubling feature of energy landscape 
theories is the inability to treat hard-sphere systems, 
where barrier crossing events are completely entropic in 
nature. In order to address this issue, we generalize the 
inherent structure concept in a manner similar to that 
discussed by Santen and Krauth 22 . Specifically, we study 
a binary hard-sphere system consisting of 225 particles 
with o\ =0.1 and and 25 particles with <ii = 0.5 at 
various packing fractions rjQ. Starting with well equili- 
brated configurations generated at a given packing frac- 
tion, the system is "crunched" until hard-sphere overlap 
occurs. Some Monte Carlo moves are made during this 
process to ensure the system reaches a stable packing 
structure. The final, stable configuration obtained from 
this procedure is called an inherent structure. In Fig. 2a 
we show the average inherent structure volume {V isijf)) 
versus packing fraction for the binary hard-sphere sys- 
tem. Remarkably, an onset packing fraction r] may be 
defined from the inherent structures that again coincides 
with the ergodic-nonergodic critical packing fraction as 
found from Eqs.(l)-(3). This not only demonstrates the 
robustness of the correlation between the dynamics of 
the ergodic-nonergodic transition as found directly from 
MCT via Eqs. (l)-(3) and the onset temperature (pack- 
ing fraction) T a (t) ), but it also hints at a deep con- 
nection between inherent structures labeled by potential 
energy, and configurations defined by packing and free 
volume. To strengthen this connection, we reconsider a 
thermal system, namely the soft-sphere system studied 
in Fig. 1. Using a crunching procedure (with an upper 
energy cutoff) similar to that used in the hard-sphere 
system, we calculate Vjs(T) versus T. In Fig. 2b we 
show that the onset temperature for changes in V/s(T) 
quantitatively coincides with T as extracted from in- 
herent structure energies (Vis(T)) for the same system. 
This demonstrates that the changes in the inherent struc- 
tures that are sampled in the liquid range where dy- 
namics become nonexponential are associated with sharp 
changes in free volume and structural packing motifs. 
In fact, such quantities are more fundamental than the 
quenched potential energy and allow for the generaliza- 
tion of landscape-like concepts to athermal systems such 
as hard-spheres. 

Despite the fact that the solution to the full set of 
Eqs.(l)-(3) yields an ergodic-nonergodic transition tem- 
perature that is very close to T a , one still expects MCT to 
give physically sensible results in the range T c < T < T Q . 
Indeed, Kob et al. have demonstrated that Eqs.(l)-(3) 
accurately account for dynamics at these temperatures 
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FIG. 2: (a) Average "inherent structure" volume vs. pack- 
ing fraction for the hard-sphere mixture as calculated by the 
"crunching" procedure described in the text. The dashed line 
indicates the location of the ergodic-nonergodic transition as 
calculated from Eqs.(l)-(3). (b) "Crunching" procedure ap- 
plied to the soft-sphere system. An energy cutoff has been 
employed to define the final IS volume. The dashed line in- 
dicates the temperature of the ergodic-nonergodic transition 
as calculated by Eqs.(l)-(3). T in this system as defined 
via ~Eis(T), T E = (L332 (see Fig. la), agrees well with that 
extracted from the VisiT) definition, TY = 0.356. 



spin-spin correlation function for a finite N system. 
Clearly, this "succeptability" measures finite N correc- 
tions to mean-field dynamics for N , N < oo. In Fig. 3b, 
we plot Xn,n'(T) versus T for several different choices 

of N and N . Remarkably, x(T) sharply increases from 
zero very close to T a . This result is not sensitive to N 
and N for a wide range of values. Indeed, this result 
demonstrates that noticeable corrections to mean-field 
(MCT) dynamics occur at T D > T c as calculated by fits 
to transport coefficients. 



(N, N') = (50, 75) 
(N, N') = (75, 100) 
(N, N') = (100, 150) 
(N, N') = (150, 200) 
(N, N') = (200, 300) 
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if the static input is calculated at a higher (effective) 
temperature 2 ^, and many of the scaling predictions of 
MCT appear to be corroborated via computer simula- 
tions in this rang o 21 ! 25 . 

The freezing of F(k, t)/S(k) at T a (and not T c ) merely 
signals the sensitivity of MCT to changes in the pack- 
ing structure of the liquid and hints at the breakdown of 
the mean-field like approximation that is inherent in the 
idealized version of mode-coupling theory. To get a bet- 
ter feeling for the nature of this breakdown, we consider 
the dynamics of the random orthogonal model (ROM) 26 . 
Specifically, we take the Hamiltonian 

H=-2^J ij a i <T h (4) 

ij 

where <7i = ±1 are Ising spin variables, and Jy is 
an N x N random symmetric orthogonal matrix with 
J a = 0. The Glauber algorithm^! is employed for dy- 
namics. For N — > oo, this model is in the same dy- 
namical universality class as the p = 3 spherical p-spin 
model. Thus, in this limit, MCT is exact. For N finite, 
1/N corrections to the MCT should appear. It has al- 
ready been shown that for finite N, the model defined 
by Eq.(4) behaves similarly to real liquids, exhibiting 
landscape properties very similar to atomistic models, 
as well as non-trivial dynamical properties such as Nagel 
scaling— In Fig. 3a a plot of Ejs{T) versus T is shown 
for a particular value of N. The generic features of this 
plot are not sensitive to N over a wide range of values 
for finite TV. An onset temperature T ss 0.72 is found 
for this model, while T c = 0.536 and T K = 0.26. In Fig. 
3b, we plot Xn n' 

(T) = f™dt\C N (t) - C N ,(t)\, where 
Cjv(i) = jf Yli=i ( CT i{t) a 'i(^)) is the disorder-averaged 



FIG. 3: (a) Average inherent structure energy vs. tem- 
perature for the ROM model. T a = 0.72. (b) 
X/v at' (T)/{Xn n' CO) versus T, where the angle brackets rep- 
resent an average over all T. Curve fitting similar to that of 
previous figures yields an approximate onset temperature of 
non-mean-field effects of T = 0.735. 

To determine which portion of the dynamics is sensi- 
tive to finite N effects, we have calculated Cjy(i) for sev- 
eral different values of N for one particular temperature 
in the range T a > T > T c . The intermediate time ("/3- 
relaxation") regime is not strongly TV-sensitive, while the 
long-time ( "a-relaxation" ) regime shows strong finite N 
effects. This observation is consistent with a recent study 
of the dynamics of the finite-sized Random Energy Model 
by Ben Arous et a&. These authors have demonstrated 
that dynamical effects beyond the mean-field limit oc- 
cur at long times, and that the longest time behavior is 
compatible with the predictions of the phenomenological 
trap modelSS,. Recent computer simulation work has also 
shown that the long-time dynamics of atomistic liquids 
is consistent with an activated-like trap model precisely 
in the range T a > T > T c . Note that in the ROM model, 
for finite N, N , Xn n' CO essentially varies continuously 
through T c but changes sharply near T Q . Meta-basin dy- 
namics in atomistic simulations also show a sharp onset 
of trap- like behavior near T , and continuous variation 
near Tc^LULlH. It should also be noted that recent work 
by Berthier and Garrahan on kinetic facilitated models 
also demonstrates that activated processes set in at T Q , 
where local dynamical heterogeneities begin to occur—. 

In this paper, we have studied a number of model sys- 
tems and demonstrated deep connections between mean- 
field theory, mode-coupling theory, and the landscape 
paradigm. The picture that emerges from this work and 
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these previous studies is that T a marks the edge boundary 
for significant barriers that result from a non-mean-field 
smearing of the dynamical transition that occurs at T c 
in mean-field systems. Barriers that would be infinite at 
T c for an infinite dimensional liquid become finite, and 
influence dynamics not just in the vicinity of T c but at 
temperatures up to T a . Transport involves activation over 
these barriers, and is trap-like at very long times. 

What, then, can one expect from the idealized ver- 
sion of MCT? The analog between 1/N corrections to 
mean-field behavior in truncated mean-field spin glasses 
and corrections to the idealized MCT suggest that 
MCT should always yield quantitative results for the (3- 
relaxation regime. The remarkable scaling predictions of 
MCT in this regime should be unaffected by the mean- 
field nature of the approximations inherent in the MCT 
approach 25 . While MCT is able to account for properties 
such as time-temperature superposition in the long-time 
a-relaxation regime, finite corrections to MCT should 
be noticeable here. These corrections are connected to 
activated (trap-like) behavior that occurs in the range 
T Q > T. Note that the Gaussian trap model also dis- 
plays time-temperature superposition in the a-relaxation 
regime. Clear signs of dynamic heterogeneities occur in 
simulations for temperatures in the range T Q > T > T c 32 . 
This localized heterogeneous motion is likely connected to 
activated, trap-like behavior. 

While we have shown in this work that there are in- 
teresting physical connections between the temperature 



T a and the breakdown of mean-field theory that have not 
been previously discussed, an open question still remains 
regarding the nature of T c as a sharp dynamical crossover 
temperature. It is instructive to note that in the p-spin 
model, T c and T a lie very near each other. Non-mean- 
field effects tend to push these temperatures apart. A 
scaled parameter r = T °^ Tc may provide a measure of 
how strong the corrections to mean-field behavior are. 
In the binary soft-sphere mixture of Barrat et al. at unit 
density r w 0.25 while for the Lennard- Jones mixture of 
Kob and Andersen studied at unit density r ~ 1. In- 
terestingly, dynamical effects such as the appearance of 
prominent "hopping" peaks in the van Hove correlation 
function G(r, t)2£ and short-time secondary maxima in 
the non-Gaussian parameter az(t) and the nonlinear suc- 
ceptability Xi{t) occur near T c for the soft-sphere mixture 
but not for the Lennard-Jones system—. Perhaps rem- 
nants of the behavior expected at T c are only noticeable 
in systems with small r values that are in some sense 
closer to idealized mean- field systems 35 . Furthermore, 
it would be interesting to determine if the value of r de- 
creases as the physical dimension of the system increases. 
A systematic study of such open questions is underway. 
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